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SIMULATION RESULTS FOR THE VITERBI 

DECODING ALGORITHM 

ByBartusH. Batson, RobertW. Moorehead, 
and S. Zafar H. Taqvi* 

Manned Spacecraft Center 

SUMMARY 


Performance predictions for convolutional decoders using the Viterbi decoding 
algorithm are presented in this report. Bit error probability is chosen as the measure 
of performance and, by using digital computer simulations, is calculated as a function 
of energy per bit per noise spectral density for various encoder and decoder param- 
eters. Coding gains based on comparisons with uncoded, coherent phase-shift-keyed 
system performance are determined for code rates of one-half and one-third and for 
constraint lengths of 4 to 8. Both hard-decision and soft-decision decoders are consid- 
ered, and bit error probability is calculated for both systematic and nonsystematic 
codes. The effect of decoder block length on decoding performance also is included to 
provide a more complete estimate of the relationship between performance and decoder 
complexity. 


INTRODUCTION 


One characteristic that has made all-digital communications links appear increas- 
ingly attractive for many applications in recent years is that error control encoding and 
decoding can be applied to achieve significant improvements in overall link performance. 
The introduction of coding into a digital link allows, for a fixed transmit (or receive) 
power level and for an allowable bit error probability, transfer of more information per 
unit time. Alternately, for a fixed information rate, the introduction of coding can pro- 
vide a reduction in the transmit or receive power level required to maintain a specified 
error probability. The exact increase in information rate that can be achieved, or the 
amount of coding gain (allowable reduction in power level) that is realizable, depends 
on the particular class of encoding and decoding technique employed and on various en- 
coder and decoder parameters that must be selected by the communications system de- 
sign engineer. 
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Convolutional codes generally are conceded to be better than block codes for many 
channels (ref. 1), particularly with respect to ease of implementation, equipment com- 
plexity, power consumption, weight, and flexibility. Various algorithms, including 
several sequential decoding algorithms (refs. 2 and 3), are available for decoding con- 
volutional codes. Sequential decoders incorporate searchback operations in the hypoth- 
esis and testing of various paths through the convolutional code tree. Thus, if a path is 
hypothesized and subsequently is determined to be in error, the decoder has the capa- 
bility to discard that path, back up, and test other possible paths. The advantage of 
such a capability is that the probability of an undetected error appearing at the decoder 
output is extremely small. However, a memory unit is required to store past data for 
possible recall. In addition, a buffer is required to store incoming bits while the de- 
coding operations (including searchbacks) are being performed. Because there is al- 
ways some finite probability that the number of operations required to decode a 
particular bit can be quite large, it is possible for the input buffer to overflow. During 
such an overflow condition, the decoder output consists either of uncorrected channel 
bits or of erasures. 

Parameters that affect the achievable error probability when using sequential de- 
coding include information rate, code rate, code constraint length, input buffer size, 
and memory size (which determines the allowable number of searchbacks) . In general, 
for fixed power levels, a smaller error probability can be obtained by decreasing infor- 
mation rate, by decreasing code rate (adding more redundancy to the transmitted se- 
quence), or by increasing any of the other parameters noted previously. The code 
constraint length used for systems employing sequential decoding is typically rather 
large (greater than 20) . 

A primary advantage of sequential decoding is that a rather large performance 
gain is achievable. However, the variable decoding time per bit that results because of 
searchbacks is a distinct disadvantage for some applications. Other disadvantages in- 
clude the requirements for a memory unit and an input buffer. Fortunately, two of 
these disadvantages (the variable decoder output rate and the requirement for an input 
buffer) become inconsequential when the data rate is low enough to allow a large speed 
advantage (computation rate per data rate) of the logic unit. Sequential decoding, there- 
fore, is a very attractive technique for use in systems having moderate data rates (be- 
low perhaps a few megabits per second) . 

The Viterbi algorithm (ref. 4) for decoding convolutional codes recently has re- 
ceived considerable attention, largely because of certain inherent advantages over the 
various sequential decoding algorithms. This algorithm has been shown (ref. 5) to be 
maximum-likelihood and, therefore, optimum for the decoding of convolutional codes. 
The primary advantage of a Viterbi algorithm decoder is speed. Whereas the perform- 
ance gain (over no coding) achievable by using sequential decoding is limited primarily 
by the information rate, the gain achievable by using Viterbi decoding is limited pri- 
marily by the constraint length of the code and is relatively independent of the informa- 
tion rate. This independence is possible because no searchbacks are required by the 
Viterbi algorithm, and only a very small logic speed advantage is required. Viterbi 
decoders have the additional advantage of operating at a fixed rate; thus, no input buffer 
is required. 

The primary disadvantage of Viterbi decoders is that the decoder hardware com- 
plexity increases exponentially with increasing code constraint length. Because 
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performance gain also increases with constraint length, hardware constraints impose 
a limit on the performance gain achievable by using Viterbi decoders. In practice, 
these hardware constraints dictate that the constraint length be limited to approxi- 
mately 8. Fortunately, the decoder hardware requirements do not increase substan- 
tially for lower rate codes or when soft decisions are provided by the receiver. 


SYMBOLS 


CD coding gain 

energy per information bit 

E c energy per channel bit 

K constraint length of convolutional code 

L search length or block length of Viterbi decoder 

N q single-sided noise spectral density 

P g bit error probability 

Q number of receiver quantization levels 

R rate of convolutional code (R = 1/V) 

Sg(t), S^(t) transmitted waveforms corresponding to binary 0 and 1 

Sg(t),S*(t) received (noisy) waveforms corresponding to binary 0 and 1 

V number of encoded (channel) bits per information bit 

FUNDAMENTALS OF DIGITAL COMMUNICATIONS SYSTEMS 
EMPLOYING CONVOLUTIONAL ENCODING 
AND DECODING 


In the following sections, the concepts involved in determining the performance 
of the general coded digital communications system are introduced, and the basic con- 
cepts of convolutional encoding and decoding are discussed. This material is intended 
merely to provide background information that may aid some readers in developing a 
more complete understanding of the convolutional decoding problem. 
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The General Coded Digital Communications System 


A simple coded digital communications system is illustrated in block diagram 
form in figure 1. The system consists of a source, which generates data in the form 
of binary digits (information bits) ; an encoder, which converts each information bit into 
V channel bits (V = 2 for this example) according to a certain scheme governed by the 

code; and a transmitter, which generates 
analog waveforms S (t) and S (t) corre- 





Sjlt) 


S](t) 




sponding to 0 and 1 for transmission through 
the channel. The signal is corrupted by 
noise (generally assumed to be additive, 
white, and Gaussian) in the channel, and the 
input to the receiver consists of a sequence 
of noisy waveforms Sj^t) and S*(t). A de- 
cision device, which processes the noise 
waveforms and provides estimates of the 
corresponding transmitted channel bits, is 
incorporated in the receiver. Because the 
input to the decision device is noisy, there 
is always some finite probability that a bit 
decision is erroneous. It has been shown 
(ref. 6) that, for binary signaling over an 
additive, white, Gaussian noise channel and 
for optimum (correlation or matched filter) 
detection, the probability of bit error at the 
receiver output is 


Figure 1.- Block diagram of coded 
digital communications system. 



(1) 


where E c is the signal energy of a channel bit at the input to the decision device, 

is the single- sided noise spectral density of the receiver, and erfc ( ) is the comple- 
mentary error function defined by 


r°° 2 

erfc X = — j e - ^ d£ (2) 

s/¥ J X 


If coding is not employed by the system of figure 1, however, the receiver bit detector 
makes decisions on information bits directly. With constant transmit and receive power 
levels for the coded and uncoded systems, more signal energy is available for a direct 
decision on an information bit than for a decision on a channel bit. The redundancy 
added to the transmitted signal when coding is incorporated into the system results in 
less energy per channel bit and, therefore, in a higher bit error probability at the re- 
ceiver output. The task of the decoder, which operates on the reconstructed sequence 
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of channel bits at the receiver output, is to correct as many bit errors as possible. 

For the achievement of a net coding gain, the information bit error probability after de- 
coding must be less than it was for the uncoded system. »An uncoded information bit 
error probability curve can be constructed by substituting (energy per information 

bit) for E c in equation (1). This uncoded bit error probability curve is shown in 

figure 2. 



Figure 2. - Bit error probability as a 

function of E,/N_ for an uncoded 
b 0 

digital communications system. 


Several encoding techniques are avail- 
able for incorporation into digital transmis- 
sion systems. These techniques generally 
are categorized as either block or convolu- 
tional. For block codes (frequently referred 
to as algebraic codes), a certain structured 
block of channel bits is assigned to each pos- 
sible group of information bits. Block codes 
are highly structured, in a mathematical 
sense, and the various decoding algorithms 
for block codes generally either exploit code 
properties that result from this mathemat- 
ical structure or apply probabilistic infor- 
mation obtained from the received signal. 

For the reasons pointed out in the introduc- 
tion to this paper, convolutional codes were 
chosen for the current investigation. 

Convolutional Encoding Fundamentals 

For convolutional encoding, each in- 
formation bit (rather than each block of in- 
formation bits) is encoded into V channel 
bits. A simple convolutional encoder, which 
consists of a shift register with K stages 
connected in some prescribed manner to V 
modulo 2 adders, is shown in figure 3. For 
each input bit shifted into the register, there 
are V encoded output bits, corresponding 
to one revolution of the commutator. There- 
fore, the length of the output (channel) se- 
quence will be V times the length of the 
input message. The rate of the code is de- 
fined as 



(3) 


and is a measure of the number of message (information) bits per transmitted symbol. 
The length of the shift register K is referred to as the constraint length of the code and 
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Figure 3. - Binary (K, V) convolutional 
encoder. 


is, roughly, a measure of the duration in 
which the encoded output bits are affected 
by any particular input bit. That is, each 
group of V channel bits depends on the cur- 
rent information bit and on the K - 1 pre- 
vious information bits. 

A particular convolutional encoder 
may be described in terms of a set of gener- 
ator coefficients that specify which stages of 
the shift register are connected to each 
modulo 2 adder. For example, two typical 
encoder configurations are shown in figure 4. 
The binary digits represent actual connec- 
tions to each adder (e. g. , binary 11101 11011 
indicates that the first adder is connected to 
the first, second, third, and fifth stages of 
the five-stage register, and that the second 
adder is connected to the first, second, 


fourth, and fifth stages). The binary generator coefficients sometimes are converted 


to octal form for notational convenience. Thus, binary 11101 11011 becomes octal 35 33, 


and binary 110101 101111 becomes octal 65 57. 



(a) First example. (b) Second example. 

Figure 4. - Typical convolutional encoder connections. 


An alternate means of describing a particular encoder configuration is the gener- 
ator sequence. For a (K, V) convolutional encoder, the generator sequence is equivalent 
to the output sequence that results from transmitting a K-bit input message consisting 
of a leading 1 followed by all zeros. One method by which the generator sequence can 
be determined is shown in figure 5. However, there is another way of determining the 
generator sequence without construction of a table. 

Note that if a single 1 is located in some stage (e. g. , stage x) of the shift register, 
the output digit from any given adder will be a 1 if that adder is connected to stage x, or 
a 0 otherwise. Thus, the first V digits of the generator sequence represent the mod- 
ulo 2 adders that are connected to the first stage of the K-stage register, the second V 
digits represent the adders that are connected to the second stage of the register, and 
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so forth. For the example shown in figure 5, 
all three adders are connected to the first 
stage; therefore, the first three digits of the 
generator sequence are 111. Because only 
the first and third adders are connected to 
the second stage of the register, the second 
three digits of the generator sequence are 
101. Likewise, the last three digits of the 
generator sequence are 100. The length of 
the generator sequence (111 101 100) is 
KV = (3)(3) = 9 bits. 

Because the generator coefficients 
specify which stages of the shift register are 
connected to each modulo 2 adder and the 
generator sequence specifies which modulo 2 
adders are connected to each stage of the 
shift register, there should be (and indeed 
is) a direct relationship between the gener- 
ator sequence and the generator coefficients. 
Given the generator coefficients, the generator sequence can be determined by inspec- 
tion. Thus, for the example shown in figure 4(a), in which the generator coefficients 
are 11101 11011, the generator sequence is 11 11 10 01 11 and is obtained by pairing 
like digits (first with first, second with second, etc.) of the two generator coefficients. 

A reverse procedure, in which the generator sequence is first separated into groups of 
V digits, can be followed to obtain the generator coefficients from the generator 
sequence. 

For those particular encoder configurations in which the first modulo 2 adder is 
connected to only the first stage of the shift register, the resulting codes are said to be 
systematic. For systematic codes, then, the first bit of each V-bit channel sequence is 
the same as the current information bit. All codes that are not systematic are said to 
be nonsystematic. 

As will be shown later, the Viterbi decoding algorithm yields better performance 
with nonsystematic convolutional codes. From a set of convolutional codes, only cer- 
tain "good” codes are used for implementation. These good codes lead to minimum de- 
coding error probability. A computer search technique usually is applied to select a 
good convolutional code, based on preselected criteria (such as maximized minimum 
distance, maximized free distance, and noncatastrophic error propagation properties). 


K-3 
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Figure 5. - Determination of the gener- 
ator sequence for a (3, 3) convolu- 
tional code. 


Convolutional Decoding Fundamentals 

Each bit shifted into the K-bit register results in one of two possible V-bit outpht 
sequences. One possible output sequence corresponds to shifting in a 0, whereas the 
other corresponds to shifting in a 1. The specific V-bit sequence that results when a 
bit is shifted into the register, however, depends on the previous K - 1 bits that are 
retained in the register. Thus, a given input message bit affects the current V-bit out- 
put sequence and the next K - 1 V-bit output sequences as well. 
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The behavior of any convolutional encoder may be illustrated diagrammatically 
by a tree structure, as shown in figure 6 for the encoder of figure 5. The labels on the 
branches indicate encoder outputs. By convention, the code tree is arranged so that the 
upper branch from any node corresponds to shifting a 0 into the K- stage register and the 
lower branch to shifting in a 1. The encoded output sequence corresponding to a given 
input message sequence may be found by following the appropriate path through the code 
tree. For example, an input message sequence of 1011. . . results in an output sequence 
of 111 101 Oil 010 


000 


Encoded 

output 

digits 

000 


000 


Input 

digits 


111 


101 


111 


010 


100 


101 


111 


Oil 


010 


001 


no 


000 


Decoding of convolutional codes in- 
variably is based on the code tree structure. 
Sequential decoding algorithms assume a 
tentative transmitted message, encode this 
message with a replica of the encoder, and 
compare the resultant coded output sequence 
with the actual received sequence. If these 
two coded sequences agree to within some 
specified amount, decoding is assumed to 
have been accomplished. If agreement be- 
tween the two sequences does not meet the 
desired criteria, another tentative message 
is assumed, and the process is repeated. 

The optimum convolutional decoder 
chooses the path through the code tree that 
has maximum likelihood, given the received 
sequence. That is, the decoded message 
sequence will provide a coded output se- 
quence that is closer (differs in fewer bit 
positions) to the actual received sequence 
than the coded output sequence correspond- 
ing to any other possible message sequence. 
Because the code tree apparently is infinite 
(i. e. , the number of branches doubles each time a bit is shifted into the encoder), 
choosing a maximum-likelihood path through a code tree such as that shown in figure 6 
would appear to be a hopeless problem. This is not the case, however, as will be 
pointed out in subsequent discussion. 
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Figure 6. - Code tree for a (3, 3) con- 
volutional encoder with generator 
sequence 111 101 100. 


Inspection of figure 6 reveals that, although the code tree does grow without bound 
as more and more message digits are shifted into the encoder, the growth is completely 
repetitive after a point. Note that as the Kth (third) message bit is shifted into the en- 
coder, eight possible output branches (000 to 111) exist. As the (K + l)th message bit 
is shifted in, there are 16 possible branches, but only eight of these are distinct. Be- 
cause the upper eight branches are identical to the lower eight branches, it is not neces- 
sary to show all 16 branches on the diagram. The tree structure could be simplified 
greatly by connecting nodes A and E, B and F, C and G, and D and H. If this is done 
for the (K + l)th input bit, it becomes evident that the same procedure could be re- 
peated for the (K + 2)th input bit, again for the (K + 3)th input bit, and so forth. Thus, 

K K- 1 

at no location in the tree is it necessary to draw more than 2 total branches or 2 
total nodes. 
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Additional insight into the basic simplicity of the convolutional tree structure can 
be gained by examining the state diagram of the encoder. An encoder state is defined 
as the contents of the first (most recent) K - 1 stages of the K-stage shift register. 

K - 1 K-l 

There are 2 possible states, corresponding to the 2 possible combinations of 

K-l binary digits. Given an encoder state, only two possible encoder states can be 
entered as a message digit is shifted into the register. The first possible state corre- 
sponds to shifting in a 0, whereas the other possibility corresponds to shifting in a 1. 
Because the most ancient (Kth) bit in the encoder register is dumped as a message bit 
is shifted in and, therefore, cannot affect any subsequent state, that Kth bit is not con- 
sidered when defining the register state. 

The state diagram for the (3, 3) convolutional encoder that was considered in pre- 
vious examples in this report is shown in figure 7. The circled numbers represent en- 
coder states, the single digits (in parentheses) represent input message digits, and the 
triple digits represent the encoded output digits that occur as the encoder state is 
changed. For example, if the current encoder state is 10, two possible states can be 
reached as a message digit is shifted in. These possible states are 01 and 11 and cor- 
respond, respectively, to 0 and 1 inputs. Assuming a 0 input, the total shift register 
contents will be 010, resulting in the encoded output digits 101. Alternately, if a 1 is 
shifted in, the register contents will be 110, resulting in the encoded output digits 010. 
This procedure was followed to obtain the complete encoder state diagram shown in 
figure 7. 




(a) Encoder. (b) State diagram. 

Figure 7. - State diagram for a (3, 3) convolutional encoder with generator 

sequence 111 101 100. 
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The state diagram allows a convenient and straightforward determination of the 
output message corresponding to a particular sequence of input digits. For example, 
if the input sequence is 10110. , and the encoder is initially at state 00 in figure 7, 
the first digit (a 1) results in an output of 111 and a new encoder state of 10. The sec- 
ond digit (a 0) changes the state to 01 and provides an output of 101. The third digit 
(another 1) changes the state back to 10 and provides an output of Oil, and so forth. 

This procedure could be repeated for an indefinite sequence of input message bits and 
is more convenient than tracing through the tree diagram. 

The state diagram can be redrawn in a form that is even more convenient for dis- 
cussing the operation of the Viterbi decoding algorithm. The redrawn state diagram is 

K- 1 

shown in figure 8, in which the 2 states are represented as nodes on a trellis 
structure. The branches of the trellis represent possible moves from state to state. 

It should be emphasized again that only two possible states can be entered from a given 
state. Likewise, a given state can be entered from one of only two possible previous 
states. 


The possible moves (from state to 
state) are independent of the particular gen- 
erator sequence being used. For example, 
state 00 can be entered only from state 00 
or state 01, regardless of the modulo 2 
adder connections of the encoder. However, 
the encoded output digits that correspond to 
each move (and are labeled beside each pos- 
sible move) are code-dependent. For in- 
stance, a different generator sequence 
results, in general, in encoded output digits 
other than 101 and 001 corresponding to the 
two paths terminating in state 01. 

In the trellis structure of figure 8, the 
uppermost of the two paths leaving any one 
node corresponds to shifting a 0 into the 
register, while the lower path corresponds 
to shifting in a 1. This statement can be 
verified by comparing figure 8 with the orig- 
inal state diagram shown in figure 7. All 
four paths terminating in the upper two 
states (00 and 01) are paths that resulted when a 0 was shifted into the register, and all 
four paths terminating in the lower two states (10 and 11) are paths that resulted when 
a 1 was shifted in because the first digit of a new register state has to be the same as 
the message digit that was just shifted in. Thus, states 00 and 01 can be entered only 
if a 0 is shifted in. 

The encoded output sequence corresponding to a given input message sequence can 
be determined entirely from the trellis structure of figure 8. However, for convenience, 
this structure is sometimes expanded to the form shown in figure 9. The expanded trel- 
lis in figure 9 actually is an alternate way of viewing the code tree shown earlier in fig- 
ure 6. Although only the first set of nodes is labeled on the expanded trellis structure, 



state 

Figure 8. - Redrawn state diagram for 
a (3, 3) convolutional encoder with 
generator sequence 111 101 100. 
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000 000 



Figure 9. - Expanded version of the 
state diagram (trellis structure) 
for a (3, 3) convolutional encoder 
with- generator sequence 
111 101 100. 


o o 


successive nodes at the same level corre- 
spond to the same state (i. e. , the uppermost 
node always corresponds to state 00, etc.). 
The manner in which a path through the ex- 
panded trellis is determined for a typical 
input message sequence (10110. . .) is illus- 
trated in figure 10. The resulting output 
sequence is 111 101 011 010 001. . . . 

The Viterbi decoding algorithm, which 
will be discussed in the next section of this 
report, is visualized best in terms of the ex- 
panded trellis structure. This algorithm 
basically attempts to find a path through the 
trellis that is as close as possible (differs 
in the fewest bit positions) to the received 
encoded sequence. The information se- 
quence (input message sequence) correspond- 
ing to this path then is assumed to be the 
same as the original input message at the 
encoder. Given the received channel se- 
quence, the algorithm is optimum in the 
sense that the most probable transmitted 
message is selected. 


DESCRIPTION OF THE VITERBI 
DECODING ALGORITHM 


'v/ From the state diagrams (including 

110 ° ° ° ^ ° the trellis structures) shown previously in 

. . this report, it can be observed that once 

two paths are in the same state (by merging 

irv . , . , ,. ,, , ,, at a common node), both paths have identical 

Figure 10. - A typical path through the , . , , .V , , ' • _ ... . 

° . , . ° extensions out of that state. Because this is 

expanded trellis structure for a ,, , . , ,, , , ,, ,, 

(3. 3) convolutional encoder with the h case - ‘ can be obse . r f' 1 that „ bolh 

generator sequence 111 101 100. subsequently will correlate equally well with 

n the received sequence. Therefore, because 

the objective is to find the path that corre- 
lates best with the received sequence, it should be possible to eliminate one of the two 
paths that enter any state. By eliminating the smaller correlated path entering each 
state, only those paths that never can be candidates for the highest correlated path 
through the trellis structure are discarded. This is the general idea on which the 
Viterbi decoding algorithm is based. 


Figure 10. - A typical path through the 
expanded trellis structure for a 
(3, 3) convolutional encoder with 
generator sequence 111 101 100. 


The Viterbi algorithm was discussed first by A. J. Viterbi (ref. 4) and was shown 
to be asymptotically optimum. In a later paper (ref. 5), G. D. Forney observed that 
Viterbi 's algorithm is optimum in the maximum -likelihood sense. 
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When implementing a Viterbi algorithm decoder, several hardware variations are 
possible. The particular implementation to be considered here maintains a running 

K- 1 

score, or branch metric, on each of the 2 most likely paths through the trellis 

K- 1 

structure. The information sequences corresponding to these 2 paths are called 
survivor sequences. At each step in the algorithm, the survivor sequences terminating 
K- 1 

in each of the 2 nodes are determined. The scoring system is such that the most 
likely path through the trellis is the path having the lowest score. 

The implementation of the Viterbi algorithm consists of the following steps. 

1. All scores and survivor sequences initially are set to 0. 

2. A received branch (V-bit segment) is correlated with each of the two possible 

K- 1 

branches out of each of the 2 states, and delta scores (AS's) are generated. A AS 
is the number of bit positions in which the received branch differs from the branch with 

which it is being correlated. There are 2 such correlations to be performed. For 
example, for the K = 3, V = 3 case that was considered previously in this report, it is 

3 

necessary to correlate each received 3-bit branch with each of the 2 =8 possible 3-bit 
branches (000 to 111). 

3. The AS's for the two paths leaving each state are added to the previous scores 
(initially 0) for that state (actually, for the path which previously terminated in that 
state) . 

K- 1 

4. Scores for the two paths terminating in each of the next 2 states are com- 
pared. The path having the lowest score is retained, and the path having the highest 
score is dropped. (In case the two scores are equal, one of the paths is dropped arbi- 

K- 1 

trarily.) Thus, 2 running scores are retained. 

K- 1 

5. The survivor sequences for the paths terminating in each of the 2 states 

then are stored, along with their running scores, and steps 2 to 5 are repeated. A set 

K- 1 K-l 

of 2 registers is required to store the 2 survivor sequences, and another set 

K- 1 

of 2 registers is required to store the running scores for each of these survivor 
sequences. A straightforward scheme for storing the survivor sequences and scores is 
to provide a survivor sequence (SS) register and a score (SC) register for each state. 
That is, the SS Q q and SC Q0 registers always will be used for storage of the survivor se- 
quence (and its score) that terminates in state 00. Similarly, the SS Q1 and SC^ regis- 
ters will be used for the sequence terminating in state 01 and so forth. This scheme 
requires that the capability exists for transfer of the entire contents of one register to 
another register. For example, if the survivor path that terminates in state 00 is the 
one from state 01, the new contents of register SSqq should be the contents previously 

stored in register SS^, plus the newest bit (a 0 in this case) that resulted in state 00 
being entered. Similarly, the new contents of register SCqq should be the contents pre- 
viously stored in register SCq^, plus the AS for the branch between states 01 and 00. 
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Operation of the Viterbi algorithm can be visualized by means of an example. The oper- 
ation for a noise-free received sequence (no bit errors), showing the encoder and its 
state diagram and illustrating the manner in which the survivor paths are traced through 
the trellis structure, is presented in figure 11. A step-by-step summary of the SS and 
SC register contents is contained in table I. Step 1 is merely an initialization step, 
during which. all the SS and SC registers are set to 0. 


3 

At step 2, the first received branch (111) is correlated with the 2 =8 possible 
branches (000 to 111) shown in figure 11(b). The AS's that are generated are equal to 
the number of bit positions in which the received branch differs from the branch under 
consideration. Thus, the AS between the first received branch and the 000 branch is 3, 
the AS between the first received branch and the 111 branch is 0, and so forth. The two 
paths terminating in each of the four states then are compared; the path with the lowest 
score is retained, and the path with the highest score is dropped. For state 00, the 

000 branch (which has a score of 3) is 
dropped in favor of the 100 branch (which 


Message 

digits 

1011000 ... 


Encoded output digits 
111 101 011 010 001 100 000 ... 



has a score of 2). Likewise, the 101 branch 



(a) Encoder. (b) Encoder state diagram. 

Figure 11. - An example of Viterbi algorithm operation (noise free). 
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Transmitted sequence 
Error sequence 


111 

000 


101 

000 


Oil 

000 


010 

001 

100 

000 

000 

000 

000 

000 


Received sequence 



111 

101 

Oil 

010 

001 

100 

000 

Step 


Step 


Step 


Step 


Step 


Step 


Step 


Step 

1 


2 


3 


4 


5 


6 


7 


8 



Figure 11. - Concluded. 


is retained for state 01, the 111 branch for state 10, and the 110 branch for state 11. 
Scores for each of the surviving paths (only one branch in length at this stage) are shown 
encircled under the states in which the paths terminate. As shown in table I, these 
scores, plus the survivor sequences, are stored in the appropriate SC and SS registers. 
(Although the scores are indicated in decimal form for convenience, they actually would 
be stored in binary form.) 

At step 3, the second received branch (101) is correlated with the eight possible 
branches to yield new AS's. The AS for each branch is added to the score for the state 
which that branch leaves, giving eight new scores. The two paths terminating in each 
state are compared, and the path with the lowest score is retained. For purposes of 
illustration, in case the two paths being compared have equal scores, the uppermost 
path is retained arbitrarily, and the lower one is rejected. Hence, of the two possible 
paths terminating in state 10, each having a score of 3, the path coming out of state 00 
is selected. Similarly, the path coming out of state 10 is selected as the path termi- 
nating in state 11. Again, the storage of survivor sequences and scores at the end of 
step 3 is shown in table I. The most ancient bits in the SS registers are those on the 
left-hand sides. The contents of the SSq Q register are obtained by shifting in the pre- 
vious contents of the SSq^ register, followed by a 0. Similarly, the SSq^ register 
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TABLE I. - REGISTER CONTENTS (SS AND SC) FOR VITERBI ALGORITHM 
EXAMPLE (NOISE FREE) 


Step no. 




Register contents 




ss oo 

sc oo 

SS 01 

SC 0! 

^lO 

o 

lT 

CO 

“n 

sc ll 

1 

- 

- 

- 

‘ 

- 

- 

-- 

- 

2 

0 

2 

0 

1 

1 

0 

1 

1 

3 

00 

2 

10 

0 

01 

3 

11 

3 

4 

100 

3 

110 

4 

101 

0 

Oil 

4 

5 

100 0 

4 

101 0 

3 

100 1 

5 

101 1 

0 

6 

100 00 

5 

101 10 

0 

101 01 

4 

101 11 

3 

7 

101 100 

0 

101 010 

5 

101 101 

3 

101 111 

4 

8 

101 100 0 
. 

0 

101 101 0 

5 

101 100 1 

3 

101 101 1 

4 


contents are obtained by shifting in the previous contents of the SS^q register, followed 

by a 0; the SS^ register contains the previous contents of the SS Q0 register, followed 

by a 1; and the SS register contains the previous contents of the SS „ register, fol- 
lowed by a 1. 

Operation of the algorithm for steps 4 to 8 of the current example is contained in 
figure 11 and in table I. The reader is encouraged to trace the various steps through 
in detail and to verify the tabulated results. The register interchange operations, in 
particular, should become much more obvious as this is done. 

Another example of the Viterbi algorithm is provided in figure 12 and in table II. 

In this case, however, it is assumed that some of the received bits are in error. Again, 
the reader is encouraged to follow the procedure that was followed in the previous ex- 
ample and to verify the tabulated results of figure 12 and table II. 

The basic operation of the Viterbi algorithm should be clear at this point. How- 
ever, one very important question has not yet been answered: When is a bit decision 
made? In tables I and n, it is shown that as each successive branch is received, the 
K- 1 

lengths of the 2 survivor sequences increase by 1 bit. It is necessary to make a 

K- 1 

decision on a bit after some finite length of time. Up to this point, all 2 survivor 
sequences simply have been stored along with their associated scores. If the SS regis- 

K- 1 

ters are made L bits long and if L is large enough, all 2 SS registers eventually 
will agree on the initial bit in the survivor sequences. For the noise-free example 
(table I), all registers agreed on the first two bit positions after only five steps. For 
the example with errors summarized in table II, however, five steps were required for 
a unanimous decision on only the first bit. As will be indicated in the subsequent simu- 
lation results, if the SS register lengths are made equal to approximately five or six 
constraint lengths (L « 5K), all registers will, with high probability, agree on the initial 
few bits of the survivor sequences, even when the data are very noisy (i. e. , there are 
many received errors). 
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TABLE H. - REGISTER CONTENTS (SS AND SC) FOR VITERBI ALGORITHM 

EXAMPLE (WITH ERRORS) 


Step no. 



Register contents 




ss oo 

SC 00 

SS 01 

SC 01 

SS 10 

SC 10 

SS 11 

sc n 

1 

-- 

-- 

-- 

— 

-- 

-- 

-- 

-- 

2 

0 

1 

0 

2 

1 

1 

1 

0 

3 

00 

2 

10 

2 

01 

3 

11 

1 

4 

000 

4 

110 

2 

101 

2 

111 

3 

5 

110 0 

3 

101 0 

4 

110 1 

4 

101 1 

3 

6 

110 00 

4 

101 10 

3 

110 01 

5 

110 11 

6 

7 

101 100 

3 

110 010 

6 

no ooi 

6 

no on 

7 

8 

101 100 0 

3 

110 001 0 

8 

101 100 1 

6 

no ooi i 

7 

9 

101 100 00 

3 

101 100 10 

8 

101 100 01 

6 

101 100 11 

7 

10 

101 100 000 

3 

101 100 010 

8 

101 100 001 

6 

ioi ioo on 

7 


Thus, a bit decision can be made after an L-bit delay in which L is of the order 
of five or six constraint lengths. To allow for the unlikely event in which all SS regis- 
ters do not agree on the oldest bit in storage, the decoding decision rule will be to 
choose the oldest bit of the survivor sequence having the lowest score. A decoded bit 
can be shifted out of one end of the L-bit register as a new bit is shifted into the decoder. 
Thus, except for an L-bit delay, decoding is performed in real time, on a bit-by-bit 
basis. 
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The oldest bits contained in the SS registers of tables I and II (after several steps) 
do indeed correspond to the original message input bits at the encoder. The second ex- 
ample (fig. 12 and table II) not only illustrates operation of the Viterbi algorithm, but 
demonstrates the error correction capability of a convolutional code that is decoded by 
that algorithm. 

An important decoder design parameter is the required survivor sequence register 
length L, sometimes referred to as the decoder search length or block length. As dis- 
cussed in the introduction, the performance gain over no coding that is achievable by 
using Viterbi algorithm decoding is dependent heavily on code constraint length K. 

Other parameters that affect the achievable performance gain, but do not impact the re- 
quired decoder hardware as severely as constraint length, are code rate R and number 
of receiver quantization levels Q. The quantization process involves representing the 
integrated signal level (from the receiver demodulator) as one of a finite number (Q) of 
possible levels. In general, two quantization schemes, known as ’’hard" and "soft” de- 
cisions, are used. The term "hard decisions” (Q = 2) implies that a threshold (usually 
0 volt) is set and that the input bit to the decoder is recognized as 0 or 1, depending on 
whether the integrated signal level is above or below the threshold. A disadvantage of 
hard-decision schemes is that all bit decisions are weighted equally, regardless of the 
relative proximity to threshold of the various integrated signal levels. This disadvan- 
tage is overcome by soft-decision schemes, which incorporate multiple thresholds. 

One common soft-decision scheme uses eight threshold levels (Q - 8). In general, 
loggQ bits are required to indicate to the decoder the relative magnitude of each chan- 
nel bit. 

In addition, Viterbi decoder operation is code -dependent. For example, it has 
been shown (ref. 7) that better performance is obtainable if nonsystematic codes are 
used rather than systematic codes. The effects of these various parameters and of 
various code types on the Viterbi decoder performance are reported in the following 
sections of this paper. 


SIMULATION DATA AND PERFORMANCE PREDICTIONS 


Parametric studies of the Viterbi decoding algorithm were performed by using 
digital computer (Univac 1108) simulations. Details of the FORTRAN IV computer pro- 
gram are outlined in appendix A, and a complete program is listed in appendix B. The 
program originally used an all- zero information bit sequence (without loss of generality) 
and selected the upper of the two paths leading to any node in the trellis diagram when 
the two scores are equal. By consistently making decisions favoring the upper paths 
(which correspond to 0 information bits), the results were somewhat biased on the opti- 
mistic side. The magnitude of this bias was determined later to be almost insignificant 
(0. 15 decibel) when the program was modified (1) to make a random path selection when 
the scores for the two paths leading into a node are equal and (2) to accept a random se- 
quence of data bits and to make random path selections. 
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Details on the various codes that were simulated are tabulated in appendix C. The 
codes used are the "good" nonsystematic codes reported in references 8 and 9 and the 
good systematic codes reported in reference 10. The data presented in appendix C in- 
clude the rate, the constraint length, the generator sequence, the generator coefficients, 
and the shift register representation of each code used in the computer simulations. 

Appendix D contains the "raw" data obtained from the Univac 1108 computer simu- 
lations. The data tabulated in appendix D for each code include the number of informa- 
tion bits used for each computer run, the value of E^/Nq (energy per information bit 

per noise spectral density) and the resulting channel error rate at the decoder input, 
the number of output (information bit) errors, and the output information bit error rate. 
All data shown on the error probability plots presented later in this section were de- 
rived from appendix D. 


Effects of Noise on Predicted Performance 

Some early simulation runs produced results that were found to be approximately 

-4 

2 decibels on the pessimistic side (at 1x10 bit error probability), when compared to 
theoretical bounds predicted by Viterbi (ref. 11). These results were obtained by using 
a random generator (RAND 4) that is based on Lehmer's method. The program is listed 
in appendix B. As indicated in the program listing, the noise is generated in two phases 
Phase I involves the generation of uniformly distributed random numbers. In phase II, 
these random numbers are translated into Gaussian noise in the dimensioned array 
QUANT. A subsequent study of various noise generation techniques disclosed that, of 
all the tests performed on the uniformity of the random numbers, the dimensionality V 
(where l/V is the code rate) of the serial and run correlation was very important for 
this application. Because simulation runs using the Ran- 10 program, which was tested 
in particular for the correlation criterion, agreed much more closely with the theoret- 
ical bounds, the Ran- 10 generator was used for the simulation data presented in this 
report. The Ran- 10 noise is generated by the power residue method (ref. 12) and has 
an initial value of 373620336005 octal and a multiplier value of 1045 octal. 


Effects of Code Rate on Decoder Performance 

The bit error probability curves shown in figure 13, which were derived directly 
from the simulation data presented in appendix D, indicate that for a constant constraint 
length, the performance improvement as a result of coding is greater for rate one-third 
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W dB 


Figure 13. - Effects of code rate on 
Viterbi decoder performance. 


than for rate one-half. For K = 5, this in- 
crease is approximately 1. 1 decibels, 
whereas for K = 8, the increase is only 
approximately 0. 4 decibel (for hard deci- 
-4 

sion) at 1 x 10 bit error probability. 

Note that for a given rate and constraint 
length, the absolute performance gain 
(over no coding) decreases with increasing 
error probability and that as the error 

probability approaches 1 x 10 \ the coding 
gain approaches 0 or even becomes 
negative. 


Effects of Receiver Quantization on 
Decoder Performance 

Figure 14 was obtained by replotting 
two of the hard-decision (Q = 2) curves of 
figure 13 and by adding a corresponding set 
of curves for soft decisions (Q = 8) . The 
soft -decision curves were obtained for the 
same set of codes as the hard-decision 
curves. Figure 14 is indicative that an ad- 
ditional performance improvement of ap- 
proximately 1. 6 to 1.9 decibels at a bit 

-4 

error probability of 1 x 10 can be realized 
by using soft rather than hard decisions at 
the decoder input. 



Figure 14.- Effects of receiver quanti- 
zation levels (Q) on Viterbi decoder 
performance. 
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V 


Effects of Code Constraint Length on Decoder Performance 

Although it can be seen from figures 13 and 14 that the performance improvement 
Obtainable by using Viterbi algorithm decoding is greater for K = 8 than for K = 5, it 
is of interest to explore this variation in somewhat greater detail. Performance curves 
for Q = 8, V = 2, and K = 4 to K - 8 are shown in figure 15(a), and the correspond- 
ing curves for V = 3 are shown in figure 15(b). These figures are indicative that the 
performance gain increases with each increase in K. As mentioned previously, how- 
ever, the hardware complexity of the decoder is highly dependent on K; consequently, 
a practical upper limit on K is approximately 8. Figures 15(a) and 15(b) are indicative 
that the increase in performance gain between K = 4 and K = 8 is approximately 

1. 0 decibel for either V = 2 or V = 3 (for soft decision) at 1 x 10~ 4 bit error 
probability. 



W® E b/ N 0- dB 


(a) Q = 8, V = 2. (b) Q = 8, V = 3. 

Figure 15. - Effects of code constraint length on Viterbi decoder performance. 
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Eb/N 0 . dB 


(c) K = 8, V = 3, Q - 8. 


(d) K = 5, V = 3, Q = 8 (systematic 
code) . 


Figure 17. - Concluded. 


CONCLUDING REMARKS 


-4 

The coding gains possible (at a bit error probability of 1 x 10 ) by using Viterbi 

algorithm decoding are summarized in table III. These coding gains were determined 
directly from figures 13 and 15 and represent gains achievable by using the nonsyste- 
matic codes. 
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TABLE III. - PERFORMANCE GAINS ACHIEVABLE BY USING VITERBI 


ALGORITHM DECODING (BEST CODES) 


K 

V 

Q 

-4 

Approximate gain achieved at P g - 1 x 10 
as compared to uncoded PSK, a dB 

Reference 

5 

2 

2 

1. 90 

\ 


8 

2 

2 

2.90 








Figure 13 

5 

3 

2 

3. 40 



8 

3 

2 

3. 70 



4 

2 

8 

3. 90 

\ 


5 

2 

8 

4. 10 



6 

2 

8 

4. 35 


> Figure 15(a) 

7 

2 

8 

4.65 



8 

2 

8 

4.95 



4 

3 

8 

4. 35 

' 


5 

3 

8 

4.80 



6 

3 

8 

5.00 


> Figure 15(b) 

7 

3 

8 

5.20 



8 

3 

8 

5.35 




a Phase-shift keying. 


The simulation results and performance predictions presented in this report can 
be used to establish certain design goals for Viterbi algorithm decoders and for com- 
munications systems employing convolutional encoding and Viterbi decoding. It was 
determined that higher performance gains (over no coding) are obtainable by using non- 
systematic codes, lower code rates, and longer code constraint lengths, and by incor- 
porating a soft-decision capability at the receiver. Furthermore, it was found that a 
decoder search length of four times the constraint length is sufficient for nonsystematic 
codes, whereas a search length of two times the constraint length was sufficient for the 
systematic code considered. 

The relatively large coding gains realizable, together with the hardware advan- 
tages discussed in this report, make the Viterbi decoding algorithm especially attrac- 
tive for incorporation into many practical communications systems. Viterbi decoding 
appears to be particularly well suited for the relatively high data rate systems charac- 
teristic of manned space flight. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, July 26, 1972 
914-50-50-17-72 
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APPENDIX A 

PROGRAM DEFINITION FOR COMPUTER SIMULATION OF THE 
VITERB1 DECODING ALGORITHM 


In this appendix, the digital computer program developed to simulate a communi- 
cations channel employing convolutional encoding and Viterbi algorithm decoding is 
described. The program is based on the requirements of the Univac 1108 computer and 
has been made sufficiently flexible to accept variations in the encoder and decoder pa- 
rameters. With the existing program capability, the bit error probability at the decoder 
output can be determined for any input signal-to-noise ratio and for any set of encoder 
and decoder parameters. 


PROGRAM INPUT REQUIREMENTS 


I nput Data 

The program has provisions for the following inputs. 

Number of passes . - The number of passes (PASS) must be supplied as data. Each 
PASS processes 324 information bits. This number is nine times the word length 
(36 bits) for the machine and is selected for output considerations. 

Number of adders . - The number of encoder modulo 2 adders (1 < V < 3) must be 
supplied as data (RATE) . 

Constraint length . - The constraint length K must be supplied as data (CL) . 

Encoder hookup connections . - Encoder hookup connections, known as generator 
coefficients (HOOKUP(l) to HOOKUP (RATE)), must be supplied as data. 

Quantization threshold levels . - The quantization threshold levels for the soft- 
decision logic at the decoder (QUANT(l) to QUANT(8)) must be supplied as data. 

Crossover probability . - The binary symmetric channel crossover probability (P) 
must be supplied for hard decisions. 
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Data Card Formats 


The data card formats are as follows. 

Data card 1 . - Data card 1 contains the following information. 

PASS: PASS is contained in columns 1 to 5. For example: PASS - 20; 
column 4=2, column 5 = 0. 

RATE: RATE is contained in columns 6 to 10. For example: RATE = 2; 
column 10 = 2. 

CL; CL is contained in columns 11 to 15. For example: CL = 8; column 15 = 8. 

Data card 2 . - Data card 2 contains the following information. 

HOOKUP(l): HOOKUP(l) is contained in columns 1 to 5 in octal. For example: 
HOOKUP(l) = 1110110 j 2 binary = 166 |g octal; column 3=1, column 4 = 6, column 5 = 6. 

HOOKUP(2): HOOKUP(2) is contained in columns 6 to 10 in octal. 

HOOKUP(3): HOOKUP(3) is contained in columns 11 to 16 in octal. 

Data card 3 . - Data card 3 contains the following information. 

QUANT(l) to QUANT(8) : QUANT(l) to QUANT(8) are contained in columns 1 to 5, 

6 to 10, 11 to 15, 16 to 20, 21 to 25, 26 to 30, 31 to 35, and 36 to 40 in decimal. For 
example: QUANT(4) = 741; column 18 = 7, column 19 = 4, column 20 = 1. 

P: P is contained in columns 41 to 45. For example: P = 1.0 percent = 0.01; 
column 43 = . , column 44 = 0, column 45 = 1. 

In addition to data cards 1 to 3, the program requires external noise subroutines 
RAN and ZOR, which generate sets of uniformly distributed random numbers "be- 
tween 0 and 1.0. 


PROGRAM OUTPUT DETAILS 
Output Data 

For the supplied data input discussed previously, the program provides the follow- 
ing" outputs. 

1. The encoder hookup connections (generator coefficients) used in generating the 
convolutional code (in octal) 

2. The quantization threshold levels corresponding to a given information bit 
signal-to-noise ratio or bit error rate 
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3. A complete printout of the quantized channel bits (three channel bits corre- 
spond to one information bit) for each PASS (Each PASS processes 324 information 
bits. ) 

4. A complete printout of the decoder output (information bits) for the all- zero 
input and for each PASS (This output is in octal; that is, if, for example, the decoded 
message is 010001000. . . , the printout will be 210.) 

5. The number of output decoding errors in each PASS and the accumulated error 

rate 

After all required bits are processed, the following information is printed out. 

1. The message length (number of information bits processed) 

2. The encoded message length 

3. The total number of decoding errors 

4. The probability of decoding errors (information bits) in percent. 

The input error rate that corresponds to the quantization threshold levels supplied also 
is printed out. 

I 

A Sample Output 

For a K = 5, V = 2 convolutional code for which the necessary parameters are 
supplied to the program as data, the sample output will be as follows. 


SOFT DECISION 3. 0 PERCENT FOR (5, 2) 

HOOKUP - 00035 00023 00000 

QUANT = 432 630 797 908 966 990 998 1000 

ERROR RATE XX 

ERRORS THIS PASS XX The total number of outputs of this 

INPUT CHANNEL BITS XX type is equal to PASS. 

DECODED INFORMATION BITS XX 

SUMMARY OF RESULTS 

Noise quantization level corresponding to 3. 0 percent (information) 


MESSAGE LENGTH 133488 

ENCODED MESSAGE LENGTH 266976 

TOTAL NUMBER OF DECODING ERRORS 672 

PROBABILITY OF DECODING ERRORS IN . 503 

PERCENT 
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FUNCTIONAL DETAILS 


Definitions and Variables 


The complete computer program, written in FORTRAN IV for the Univac 1108 
machine, is listed in appendix B. Comment cards are inserted at appropriate places 
to make the program self-explanatory. However, the following functional definitions 
and notations are provided for the convenience of the reader. 


PASS. - PASS is the integer value of multiples of 324 information bits to be 
processed. 

RATE . - RATE is the number of modulo 2 adders of the convolutional encoder. 
CL. - CL is the constraint length of the code. 

HOOKUP. - HOOKUP represents connections (octal representation) of modulo 2 
adders to the stages of the convolutional encoder shift register (fig. A-l). 


Generator sequence ■ 111 101 110 



HOOKUP! 11 ■ no 

HOOKUP! 2) ■ 101 Generator coefficients ■ 110 101 111 

H00KUP13) - 111 

Figure A-l. - Convolutional encoder 
example. 


QUANT . - QUANT represents the 
quantization threshold levels corresponding 
to a given bit error rate. 

REG, SREG. - REG and SREG are 
equivalent variable names that represent a 

g 

maximum of 1 x 2 dimensions or a maxi- 
4 

mum of 1 x 2 storage dimensions used at a 
time of processing a single information bit. 

g 

The dimension 1x2 corresponds to a value 
of CL = 8 and can be changed accordingly. 

4 

The two 1x2 dimensions are used alterna- 
tively for two consecutive information bits. 
REG is a 36-bit sequence of 0 and 1 storing 
the decoded information bits for a block 
length of 36 bits. For a given CL (= K), the 
dimension corresponds to twice the number 
of nodes in the trellis diagram. 


8 

SCORE . - SCORE has a maximum of 1x2 dimensions. As in the case of REG, 
a maximum of half its dimension is used alternately for two consecutive information 
bits in calculating the score in the transition from one state to another (moving to a 

4 

node in the trellis diagram) . A maximum of 1 x 2 dimensions (corresponding to the 
total number of nodes) is used to store the updated score for the corresponding paths 
for deciding on the best estimate of the input bit. 
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ENCODE . - ENCODE is the channel bit input to the decoder. This input simulates 
the noise added to the channel bits corresponding to all-zero information bits. ENCODE 
has a dimension of 6 by 324, half of which is used for each alternate PASS. 

TABLE. - For a given hookup connection, TABLE represents the output of the en- 
coder for a state transition. The TABLE values for the example are: TABLE(l) to 
TABLE(8) = 000, Oil, 101, 110, 111, 100, 010, and 001, respectively. The first 
K- 1 

2 dimension of the TABLE refers to the outputs corresponding to a 0 information 

K- 1 K 

bit input, whereas from 2 + 1 to 2 correspond to a 1 information bit input. To 

the previously stored scores, the difference of the received segment and the TABLE 
values leading to a node are added, and the lower of the two scores and the correspond- 
ing path are stored. 

AMASK(l) to AMASK(36) . - AMASK is a 36-bit integer value such that 

AMASK(l) = 100000000000000000000000000000000000 


AMASK(12) = 000000000001000000000000000000000000 


The dimension refers to the position of 1 in the all- zero sequence. 

OMASK(l) toOMASK(36) . - OMASK is a 36-bit integer value such that 

OMASK(l) = 011111111111111111111111111111111111 


OMASK(12) = 111111111110111111111111111111111111 


The dimension refers to the position of 0 in the all-one sequence. 

BUFFER . - BUFFER is a 36-bit decoded message in octal representation. 

OUTPUT . - OUTPUT is the decoded information bit corresponding to the path with 
the smallest score. 

FIRST, SEC. - The FIRST and SEC pointers are set alternately to 0 and 128 (cor- 
_ 

responding to the maximum 2 ). While processing odd-numbered information bits, 
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FIRST is set to 0 and SEC is set to 128; whereas for even information bits, the settings 
are reversed. This interchange of FIRST and SEC values helps in keeping the SREG 
and SCORE values. 

TOP, BOT . - The TOP and BOT pointers are set alternately to 0 and 324. While 
processing odd-numbered blocks of 324 bits (1, 3, 5. . . , corresponding to running PASS 
values of 36, 34, 32. . .), TOP is set to 0 and BOT is set to 324; whereas for the even- 
numbered blocks, the settings are reversed. 

TOPP, BOTP . - The TOPP and BOTP pointers are set alternately to 0 and 324 
after every output printout corresponding to a PASS (324 information bits). 

INDEX . - INDEX is the counter for the output. Its value is incremented by 1 for 
each information bit processed after the first 36 bits. The output is printed out when 
INDEX = 324, then INDEX is reset to 0. This means that there is no output until 
324 bits (one PASS) plus 36 bits are processed, after which the first 324 bits are printed 
out. Subsequently, INDEX gives 324 bits output for every 324 bits processed. 

ERROR, ACCUM . - The ERROR counter counts the number of errors in each 
PASS, and ACCUM is the accumulated updated number of errors. 

A. - The A is the decoded information bit error rate in percent. 

SR, LB. - The SR and LB pointers are used to set the information bits in the 
SREG and also to calculate the error corresponding to the path with minimum score. 

H. - The H is the node pointer for the path with the smallest score. 


Block Definitions 

Although the program is not grouped into definite, discrete sections, it may be 
broken into the following blocks for the purpose of functional details. 

Block A . - Block A defines the dimensional statement and reads and writes the 
data input. 

Block B. - Block B defines the constants and calculates the OMASK and the 
AM ASK. 

K CL 

Block C . - Block C calculates the 2 (=2 ) TABLE values, which are the en- 

coder output for all possible state transitions. 

Block D . - Block D generates random numbers uniformly distributed between 
0 and 1. 0 by calling an external subroutine RAND. If soft decisions are used, block D 
generates values of ENCODE, the soft-decision channel bits corresponding to a required 
bit error rate. 

Block E. - Block E processes 324 bits for each PASS. The processing involves 
K- 1 

comparing the 2 pairs of TABLE values with the corresponding ENCODE values, 


32 



updating the REG and SCORE values for the survival paths, calculating the decoding 
error, and packing the decoded output in BUFFER. The functions of five subblocks in 
block E can be explained as follows. 

Block E^: Block E^ compares the two paths leading to a node (e. g. , TABLE(l) 
and TABLE(2) leading to node 00) with the ENCODE and updates the SCORE values. 


Block E^ 


Block E 2 updates the INDEX values for printout. 


Block E^: Block E^ packs the de- 
coded output into BUFFER and calculates 
and prints out the updated error rate, the 
number of errors, the ENCODE values, 
and the BUFFER for every 324 bits 
processed. 

Block E : Block E^ sets the SR and 
LB pointers. 

Block E,.: Block Eg updates the 

SREG and SCORE values for each bit 
processed. 

Block Eg: Block Eg resets FIRST 

and SEC values after each bit is processed 
and sets pointer H for the path with mini- 
mum score to be used in block Eg. 

Block F . - Block F resets the TOP 
and BOT pointers after each PASS, checks 
that the bits corresponding to every PASS 
have been processed, computes the final 
error rate, and prints out the summary of 
the simulation results. 

Figure A- 2 is a flow chart for com- 
puter simulation of the Viterbi decoding 
algorithm . 



(a) Start. 

Figure A- 2. - Flow chart for computer 
simulation of the Viterbi decoding 
algorithm. 











Block 

A 




(c) End. 


(b) Continuation. 


Figure A-2. - Concluded. 











APPENDIX B 

LISTING OF VITERBI DECODER SIMULATION PROGRAM 


i • 
2* 
3* 
s » 
s* 

4* 

7. 

8* 

9. 

10* 
1 1 • 
l 2 * 
13* 
19* 
IS* 
J 4 • 
17* 

1 8* 
IV* 
20* 

2 1 • 
22 * 
23* 
2 V * 
2S* 

2t* 
27* 
28* 
2 9* 
30* 
31* 
32* 
33* 
39* 
3S* 
34* 
37* 
38* 
39* 
90* 
9 | » 
92* 
9J. 
99* 


C VITERBI DtCODER ALGORITHM COMPUTER PROGRAM LISTING ON UNIVAC 1108 

C 

C BEGINNING Of BLOCK A, OEf[NES DIMENSIONAL STATEMENTS 

C HEADS AND 4R1TE5 ThE DATA INPUT 

implicit inte&eria-i) 

REAL SPRED.fACC.fiNlT.fPASS.A.AK.RANDA 

heal PDCE.P 

real rnum . random i 2oR 

ABNORMAL RANDOM, 20R 
DIMENSION MSGI 18l .OUTtRRIVI 

DIMENSION A M A S K I J 4 } , 0 M AS A ( 3 6 I , G U A N T ( B I , R E G ( 2 S 4 , 2 I 

1 . OUTPUT (498 I ,SC0R£(2S4I iBUffER(9» ,TAOleI2S4) , HOOKUP! I 0 > , 

2 ENCODE ( S8 9 3 ) , SREG ( 2S4 , 2 I 

equivalence ire&iSre&i 

data OMASK ( l 1 / 0900000000000/ • AMASM I • / 03? 7 7 77777 77 7/ 

DATA MSTART/0lO9OS/iSSTAKT/Oj73420334Qasf 
cthu program is machine independent eacept fo k «ord length changes 
cahich do affect the simulated shift register lengths and formation 
cdaTa i nput 

c ..••••••<<..< 

c 

READIS.430 I InIT.RaTE.CL 

READ ( S . 43 1 I I HOOKUP I I I , I « 1 , R a TE I 

4RlTE(4,49QllMOO*Upi|), I*|, RATE) 

READIS, 430141 ,*2 

I READIS.430.END»7SO|(«UANT(||,I.i,e|,P 

*RlTEt4,4HlHSUANlU|,l»l,6) 

PASS * INIT 
*3*4 1 *42 

PRINT 4S0, 43,41,42 

4 S 0 FORMAT II HO , 7 X , ISMPaTmS OF L E N& T M , | 3 ■ 2 S H B|TS ARE 6 E ] N G USED 4ITM, 

I 1 3 , I 8 H bits IN *URO 1 AND , 13,1 SH BJTS In 40RD2.) 

C * * 

CCOnsTANT DEFINITION 

c * * * 

C BEGINNING OF BLOCK B, CALCULATES the 

C CONSTANTS, OmaSk And AMASK 

accum.q 
f |RST*U 
SEC ■ 1 28 
TOP«Q 
801.329 
TOPP.Q 
BOTp.329 


MS* 
9 6 * 
M 7 • 
96* 
99* 
SO* 
Si* 
S 2 • 
S3* 
E 9 • 
SS* 
S 4 • 
S 7 • 
S8* 
S9* 
40* 
4 1 . 
42* 
43* 


ERR70T« 0 
I N 0 E A • U 
OMASK I 34 I • l 
DO 10 J« 1 . 39 
K ■ 3 4 - J 

10 OMASK (K. I » 2*0MASK ( K* l | 

DO 30 J ■ 2 ,34 

30 AMASK(J)«0R((AMaSK(J-1I'0MASMJII,0MABK(U-1}) 

&MASK.7 


c .............. 

CCALcULA T I ON of control constants 

c 


PA r HS«2. • ( cc- I I 


c 

C ENO OF BLOCK B 

c * 

CBUIlOING OF NODE PaTh TABLE 
C BEGINNING OF BLOCK c, calculates the 

c table values for state transitions 



69 ® 
65 ® 
66 ® 
if 
68 ® 
69 ® 
70 ® 
7l® 
72 ® 
7 3 * 
79 ® 
75 ® 
76 ® 
77 ® 

7 8 * 
7 9 ® 
80 * 
81 * 
82 * 
83 * 

6 9 ® 

85 ® 
86 * 
87 * 
88 * 
89 * 
9 0 ® 
9 1 • 
92 * 
9 3 * 
99 ® 
95 ® 
96 * 
9 7 ° 
98 ® 
9 9 ® 

1 ua® 

I 0 | • 
102 ® 
103 ® 
I 0-*® 
105 ® 
1 06 ® 
107 ® 


I 10. 
Ill® 
I 12® 
l I 3® 
1 1 9® 
I IS® 
I 16® 
117® 
l 18® 
119® 
I 20® 
121 ® 
122 ® 
123® 
12 9® 
l 25® 
126 * 
127® 
1 28® 
I 29® 
1 30® 
131® 
I 32® 
133® 
1 39® 
l 35® 
I 36® 


l • • 

KK® 2 ®P*THS 
00 90 J® I • KK 
M -0 

J J ■ * 85 ( J- I I 
00 80 N=1 iRATE 
K=AN0< JJ.MOOKuPlNI I 
KOUNT-O 
00 50 L®1 .36 
KAOO-ANOI* ,0 Ma5k IE) ) 

IF ( K A 0 D 1 50 , 50,90 

90 KOUNT*X 0 R ( KOUnT , 1 ) 

50 CONTINUE 

K®KOUNT 
L L ® 3 6 
M • M * 8 

00 70 MM® 1 , 3 
IF ( K 1 70.70,60 

60 M • OR ( M , OMASK l LL 1 | 

70 LL-LL -1 

80 CONTINUE 

90 TABLE(J)«M 

C 

C ENO OF SLOCK C 

C 

c 

C RANDOM bit generation AND QUANTIZATION 

c 

C BEGINNING OF BLOCK 0 , T h 1 5 BLOCK IS TU 

C BE CHAN&EO amen changing the 

C SIMULATION from hard decision to soft 

c decision or vice-versa, block generates the input 

C SEQUENCE to the OEcoOER AHiCh IS TmE EncUDED message 

c corresponding to random input plus noise 

DO 100 J®l .PATHS 
REG I j, l 1-0 
REG ( J , 2 I ®0 
REG ( J® I 28 . 1 I «0 
REG ( j® I 28 , 2 1 ®0 
SCORE I J I ®0 

100 SCORE! j« 1281-0 

SR"* l ♦ 1 
S*»l 
SX -2 
L** 1 

I S’SSTART 
M I .MST ART 
5 m I F f * 0 

RANOK = RanoOM 1 309^8 I 25 I 
I S 3 - 2*2139 
IMA • 282729 
110 CONTINUE 
IN- 1 
I M • 1 
I J »0 
TP ®0 

1 F l TOP . NE . 0 I TP *9 
DO 131 1 * 1,329 

RNUM = 2 0 R ! 0 ) 

IF (RNUM-.S I 1 35 * 1 36 , 1 36 

135 MSG! IM*TPI®aNo!MSG( IM®TP 1 ,AMASk( INI 1 
RMSG • 0 

GO TO 137 

136 MS G ( | M ♦ T P I . OR(MSG( IM*TP) .OMASKI IN) I 
RMSG ~ 1 

137 I N • 1 N ♦ I 

IF ( IN.LT . 37 ) GO TU 138 
IN • | 

I M • I M ♦ 1 

138 SHIFT • SHIFT /2 ♦ R M SG ® 2 • • I C L - I I 
00 130 3 « 1 ,RATE 

K®AN 0 (SHIF T , M 0 0 K U P ( J 1 I 
KOUNT * 0 
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137* 


00 8 S 0 L * 1 1 3 6 

1 i8* 


KADD • AND ( K , OMASK ( L 1 1 

1 3 9* 


IF (KADOI8SO.8SO.8Mo 

1 ‘♦0* 

8*40 

KQUNT • AORIKOONT.il 

1 M 1 ♦ 

850 

CONTINUE 

1 H 2 • 


RMSG • KOUNT • 7 

it 1* 


1 J * 1 J ♦ 1 

l **«*• 


CALL RAN ( I S , M 1 , AK I 

IMS* 


K ■ AK • 1000. 

1 **6* 


00 120 L* 1 1 8 

117* 


M • K-UUANT(L) 

IMS* 


I F (Ml 130, 120, |20 

1 H9* 

12 a 

CONTI NUE 

ISO* 


IF (L.GT • 8 1 L ~ 6 

1S1 • 

1 30 

ENCODE ( RATE*TOP* 1 0 I • AOR 1 ABS ( L- l 1 , RMS(, , 

152* 

131 

CONTI NUE 

1S3* 

C 


1 SM* 

c 

end of block d 

ISS* 

c 


l Si* 



1 S 7 • 

cinpuT to viterbi al&okitmii 

i be* 

c 


159* 

c 

beginning of BLOCK E, THIS Block PROCESSES 

160* 

c 

32M INFORMATION bits corresponding TO Each 

I6|* 


S P H E D * 0 

1*2* 


ehror-u 

1 S3* 


00 mho nnn* i , 32m 

1 6S* 


MMM.RATE*(NNN-1 1 

16S* 

c 


1 6 6 * 

CVITeRBI algorithm 

167* 

CNNN 

= NOOE COUNTER 

166* 

CJ - 

path counter i i-paThsi 

1 6 9* 

CK = 

general VARIABLE 

170* 

CL - 

ORIGINAL state COUNTER 

17 1* 

CM = 

double ORIGINAL CUUNTER 

172* 

CN = 

FINAL STATE COUNTER 

l ? 3* 

CJJ = 

holder of possible Encoder 

1 7M* 

CMC 

2 ROKKING value 

1 7S* 

CLL 

= SCORE HOLDER - FIRST 

176* 

CC MM r SCORE HOLDER - SECOND 

177* 

C 


178* 


DO MOO J* 1 , PA THS 

179* 


L = J-| 

1 80* 


M=2*L 

181* 


M = M ♦ 1 

1 8 2* 


N“ H 

183* 


NN=N-PATHS 

18m* 


IF (nN) IS0.IM0.IM0 

ies* 

IHQ 

N-N-PATHS 

186* 

iso 

JJ*T ABLE ( M ) 

187* 

c 

BEGINNING OF SUBBLOCK El. COMPARES THE T*0 

188* 

c 

PATHS LEADING to a node (KITH THE INPUT 

189* 

c 

from BLOCK 0 and UPDATES the score 

190* 

CFIRsT comparison 

19 1* 


LL-0 

19 2* 


NDK z RATE*TOP*MMM*RATE*l 

193* 


DO ISO N N ■ I , R A T E 

19m* 


KKzANDl JJ.GMAS* 1 

1 95* 


KK z KK-ENCODE ( NOA-NN 1 

194* 


IF ( KK 1 160,170,170 

197* 

160 

KKz.kk 

198* 

170 

ll=ll*kk 

1 99* 

180 

JJ; JJ/B 

2U0 • 


JO = T ABLE (M* l 1 

2 0 | • 

CSECnNO COMPARISON 

202* 


MM = 0 

203* 


00 210 nn * 1 .rate 

'2 0*4* 


kk=and( jj.gmaski 

2QS* 


KK r KK-ENCODE ( NOA.NN 1 

2 06 * 


IF ( KK I 190, 200 , 200 

207* 

1 90 

K K - ■ K K 

208* 

200 

MM=MM*KK 

209* 

310 

J J« J j/ 8 
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2 1 00 


LL*LL°SCOHE (N.F 1 R 3 T 1 

211* 


MM«MM.SCOR£IN.t.F 1RST ) 

2 i 2o 

C S T or A G£ OF OUTPUT a 1 T 

2 1 3° 


IF III 220,220.390 

2 1 so 

220 

NN* 1 N t T f PASS 

2 1 bo 


If (NN1 230.230.29o 

2 1 6° 

230 

NNoNNN-63 

21 7o 

c 

BEGINNING OF SOBBLOCA £2, UPDATES THE 

2 18° 

c 

INDEX VALUE FQK THE PKJNTqUT, STAkTS 

2 19° 

c 

COUNTING ONLY AFTER *3 BITS ARE PROCE5SE0 

2 2 0 s 


IF ( nn) 300,3oO,^9o 

22 1 o 

2 HO 

A = A N0IREG(H,L*I , OM a Sa < Lb 1 1 

222 o 


1 N u E X o INDEX. 1 

22 3° 


OUTPUT ( lNUEX-TOPPl.A 

2290 


NN° 1 NDEX'329 

2 25 0 


IF (NN) 300,250,230 

226» 

2&o 

l NDt g o q 

22lo 

c 


22 8o 

c 

END OF SuBBlOC^ 52 

229o 

c 


2 30 o 

C FORMATTING FOR PRINTOUT 

23 10 

c 


2.12 0. 

c 


2330 

C 1 J 

° first output to process 

2 3 M o 

C 1 A 

» last output bit to process 

23So 

C I N 

° increment position inoex 

2360 

c 


237° 

c 

BEGINNING OF sUaBEOCA E3, PACAS TmE DECODED 

2380 

c 

data into buffer, calculates and prints 

23 90 

c 

OUT UPOATEO ERROR rate, no* of errors 

2900 

c 


29 1 o 

c 

encode ano buffer 

292“ 


I JO 1 

2930 


1 A = 3 6 

2HM o 


TP 3 o 

29So 


IFIB0T.nE.0ITP»9 

2960 


00 290 IM-1 ,9 

2H 7o 


1 No 1 

2*400 


DU 280 I « I J , 1 A 

2H9o 


IF (OUTPUT! I ♦ T 0 P P I 1260,270,260 

2bO» 

260 

BUFFER! |MI°OR(BuffER( I M 1 , OMASA ( INI 1 

2510 


IF (FLO KIN- 1), 1,MSg(IM*TP|!.NE. lltRROR. ERROR ♦! 

2520 


GO TO 280 

253o 

270 

BUFFERl IMInA NO (BUFFER! [Ml , AMASA l l N 1 | 

2bH o 


IF (FLO! ( IN-i 1 , l ,R3 g( IM.TP) I ,NE,OIERRoR*tRRoR*l 

255“ 

280 

I N - I N o 1 

2560 


I J — I J o 3 6 

2570 


OUTERR! 1 M > o xOR(BUFFERIlM) ,MSGI JM.TPI I 

2580 

290 

1 K= 1 K« 36 

2 & 9 o 

CPR 1 nTOUTP 

2600 

600 

ACCUROACCUM.ERHOR 

26 10 


F A t C = 6 C C VJM 

2620 


F I N l T= 1 N IT 

263o 


fpass-pass 

2 6 *4 o 


a°facc/( (FinIt-fpassi ° 329 . di 

2 6 & 0 

56 1 

F0RMATI20H TRANSMlTTLO MESSAGE) 

2660 

55 1 

format i 1 gh oecooEo messaGci 

267o 

571 

formatusm error In decoded mlssagei 

2680 

c 


269° 

c 

END OF SUBBLOCK L3 

2700 

c 


27 1» 

CSHlfT OF OUTPUT REGISTER 

272° 


IF ( TOPP 1291,291. 292 

27 3° 

29 1 

T0PPo329 

2 7 H o 


80TP°0 

2 7 b o 


GO TO 293 

2 760 

29 2 

TOPPoQ 



2 77* 


B0TP»32M 




278* 

293 

CUNT iNUE 




2 7 9 • 

C 

BEGINNING OF SU88L0CK £.9, SETS THE 




2 B 0 • 

c 

SR AND Lb POINTER*, Lb FOLLOWS SR 




28 1* 

30C 

SR-SR- 1 




2 B 2 * 


IF (SR) 310,3)0,320 




2 b 3 • 

310 

1 F ( So , NE . 1 1 00 TO 3 j 1 * USE 

ONLY 

ONE. 

WORD 

2 B M * 


SK*« 1 




2 B b * 


S A * 2 




28 6* 


Sn ■ 1 




2 B 7 • 


SO TO 320 




288* 

3 1 1 

SK-&2 




2 B 9 * 


S** 1 




290* 


S A«2 




29 l • 

320 

LB-Sk- 1 




if 2* 


IF (LB) 330,330,390 




293* 

330 

IF (La. NE. 1 ISO TO 331 » USE 

ONLY 

ONE 

WORD 

299* 


L«* l 




29S* 


LB*A 1 




298* 


GO TO 3 M 0 




2 9 7* 

33 1 

L«*2 




298* 


L B ■ w 2 




299* 

c 

END OF SUbbLOt* LM 




3 UQ * 

c 

SCORE COMPARISON 




30) • 

c 

BEGINNING OF SUBBLoCK es, for E V t r t 




3 U 2 • 

c 

bit prucesseo, uPOaTES 2**K VALUES 




3 0 3* 

c 

OF SREG AN0 SCORE 




3 0M* 

3 M 0 

NN*LL-MM 




3 OS • 


1 F ( Nn 1 3 S 0 , 3 b 1 ,360 




306* 

3b 1 

CALL R A N 3 1 1 S3 , 1 M3 , AK | 




307* 


1FIAK , LT • ■ S 1 GO To 3 6 U 




308* 

C F 1 R s T CHOICE lb SMALLER 




309* 

3&0 

SR£G(J*SEC,S*)*ANU(SKtG(N. FIRST, BA) , A M A S M B R I 1 


3 1 0* 


SR£G(J*S£L,SX)*SREgIN*F1RsT,SA| 




311* 


SCORE ( J*SEC »*ll 




3 1 2* 


GO TO 370 




3 13* 

CLAST 

CHU1CE IS SMALLER 




3 1 9* 

360 

SREGtJ»SEC,S*)*ANO(SREG(N.l*f IRST , Sa I , 

AM AbK ISR) ) 


315* 


SREgiJ*SEC,Sxj*SREg(n*1*F IKST ,SA) 




3 16* 


SCORE ( J ♦ S E C 1 *MM 




3 17* 

370 

NN*m-n 




3 18* 


IF (NNI 3B0, 390,380 




3 19* 

380 

SREG(J*SEC,S«I*0R(SREG1J*SEC,S>»I , 0 M A S K I S N 1 1 



320* 

390 

CONT 1 NOE 




32 1* 

C 

END OF SUBBLOCN ES 




322* 

C 

ENU OF LOOP Of STATES 




32 3* 

MOO 

CUNT 1 NUE 




3 2 M * 

c 

INTERCHANGE OF REGISTERS 




325* 

c 

BEGINNING OF SObBLOCK E6, RESlTS THE 




326* 

c 

'FIRST' AND 'SEC' VALUES after EvERT bn IS 



327* 

c 

PROCESSED. SETS The POINTER M 




328* 

c 

CORRESPONUING TO Th£ MJnImUM SCORE 




329* 

c 

PATH 




330* 


IF (FIRST) 920.910,920 




331* 

M 1 0 

F IRST* 1 20 




332* 


SE C *0 




333* 


GO TO M 3 0 




3 3 M * 

M 2 0 

FIRST-0 




3 3 S • 


SEC. I 28 




336* 

930 

CONT 1 NOE 




337* 

CF 1N0 

M 1 N I MUM 




338* 


F *999999 




339* 


DO M 6 0 G« J .PA T *b 




3 M 0 * 


NN*e-ScORl < &*F 1 « iT 1 




3 M l * 


IF (NN| 960,960,980 




392* 

9S0 

F *SC0HE I G*F 1RS> 1 




3 M 3 • 


H«G*F 1 RST 




3 M M • 

M 6 0 

CONTI NOE 




3 9 & • 

C 

ENO OF SOBBLOC* E 6 AND BLOCK E 



* * * * 







3 M 7 * 

c 

nomerical spreao of scores 



* * * 







J M9 • 

c 

ENO OF LOOP OF 329 BITS PROCESSED 






3S0* 
3SI* 
3S2* 
3S3* 
3b*)* 
3SS • 
3S6* 
3 S 7 • 
3 b 8 • 
3S9* 
3 6 0* 
3 6|* 
362* 
363* 
3 6 9* 
36b* 
36 6* 
367* 
368* 
369* 
370* 
371* 
372* 
373* 
3 7 9* 
37b* 
376* 
3 7 7 * 
3 7 8 * 
379* 
380* 

38 1* 
382* 
383* 
389* 
38b* 
386* 
387* 
388* 
389* 
390* 
39( • 
392* 
393* 
3 9 *4 • 
39b* 
396* 

39 7* 
398* 
399* 
9 U0 • 
9 U 1 • 
9 U 2 • 
9 U 3 • 
90S* 
9 U b • 
906* 
9 U 7 * 


I • 
2* 
3* 
9* 
b* 
6* 
7* 
B* 


9 *4 0 CONTINUE 

C 8 E 0 1 N N 1 N 8 OF SLOCK f , RESETS TOP 

C AND SOT POINTERS FoR EVENT PASS, PRINts OUT 

c the final suhmart 

C ENO OF BLOCK F 

C 

IF (TOP) 970,970,980 
970 CONTINUE 
T OP* 329 
8 0 T «0 
80 TO 990 
980 CONTINUE 
TOP«0 
801*329 

990 CONTINUE 

00 bOO 1*1 .PATHS 

SCO SCONEll*FlRST)*SCOpEll*FlRST|-F 

c • • * * * • • • 

CRETnRN FOR ANOTHER PASS 

C 

ERNTOT*ERRTOT. error 
PASS* PASS* 1 
IF (PASS) b 1 0 , b 1 0 , 1 1 0 
bio AN 1 TE 4 6 , 700 I 

700 FORMAT 4 1 H 1 , // / I l 6A , 1 BHSUHmANT OF RESULTS) 
aRITE(6,701) P 

701 FORMAT 4 1 MO , 20* I 92H N 0 1 S£ QUANTIZATION LEvEl CORRESPONDING TO ,Fr.|, 
1 9 M PERCENT.) 

NMES08 •( INI T-l 1*329 
NCHNOG . RATE.NMESDG 
*R I TE 4 6 , 702 I NmESOG 

702 FORMAT 4 1 MO , 20* , I 6HMESSA8E LENGTH • . I 3 V ) 

AR I TE 4 6 , 70 3 I NCMNDG 

703 FORMAT 4 lHO, 20*,29H t NCOOEO MESSAGE LENGTH *,131) 

ARITEtA.TOAIERRTOT 

706 format 4 1 MO , 20* .33MT0T Al NyMBtR OF DECODING ERRORS .,122) 

POCE * FLOAT 4 1 00*tRRT OT ) /F LOAT 4 NMEbDG ) 

AR I TE 4 6 , 70 7 I POCE 

707 FORMAT 4 1 HO , 20* , 93MPR0BAB IL I T Y OF DECOOjnG ERROR, )N PERCENT ». 
1FI2.3I 

GO TO l 

7 b 0 C0NT1N0E 

PRINT 692.IS 

692 FORMAT I 1M0 ,b* , 39MR E INI T I AL I ZAT I ON VALUE FOR RAN10 *,0121 


c 

STOP 



c 

b 2 0 

FORMAT (ISH ERROR rate • 

,E lb.81 


b3Q 

FORMAT (2IH ERRORS this 

P A 5 b * , 

lb) 

b M 0 

FORMAT (2*,b»9A,Z0llll 



b50 

FORMAT (2A.bt9A.U20) 1 



6 1 0 

FORmaTIBM PASS • tl9,s*,6h CL * , 

I 3 , b A , 0 H RATE 

6 | 1 

FORMAT 4 2SH TOTAL DECODING 

errors 

• ,16) 

6 30 

FORMAT (8Ib.Fb.3l 



63 ) 

FORMAT! 80b 1 



6 HO 

FORMAT 4 1 MO , 7A , 9HH00KUP * 

,8IUb,2A) ) 

6 H 1 

FORMA T ( 1 HO ■ 7 * , 8 hQU an T " , 

8 4 I b I 1 



END 


SU8R0UT I Nt RAN ( 1 S , I M , A I 
1 S« I M* I S 

I IMS/ 393S97 3836 7 

IS.IS-I 11*39359738367) 

A*ABS( IS) 

A*A/393S9738367 

RETURN 

END 



C BEGINNING OF ALTERNATE VERSION OF BLOCK D. 

C THIS VERSION, WHICH IS TO BE USED FOR 

C HARD-DECISION SIMULATIONS, SHOULD 

C REPLACE STATEMENTS 115 TO 154. 

110 CONTINUE 

NN=RATE*324 
DO 130 J= 1 , NN 
CALL RAN(IS, MI, AK) 

IF(AK. LT. P) GO TO 133 
ENCODE(RATE*TOP+J)=0 
GO TO 130 

133 ENCODE(RATE*TOP+J) = l 

130 CONTINUE 

C END OF ALTERNATE BLOCK D 



APPENDIX C 

DETAILS OF CONVOLUTIONAL CODES USED IN VITERBI 
DECODER SIMULATIONS 


This appendix contains configuration details of "good” convolutional codes that 
have been suggested by several investigators (refs. 8 to 10). The parameters (K, V, 
generator sequence, and generator coefficients) of each of the various codes are sum- 
marized in tables C-I to C-III, and the same information is presented in pictorial form 
(shift register representation) in figures C-l to C-3, respectively. Although simulations 
were not performed using all the codes described in this appendix, details on the unused 
codes are included for reference and potential use by some readers. 


REFERENCE 


C-l. Heller, J. A.: Sequential Decoding, Short Constraint Length Convolutional 

Codes. Jet Propulsion Laboratory Space Program Summary 37-54, vol. Ill, 
Dec. 1968, pp. 171-177. 
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TABLE C-I. - NONSYSTEMATIC CODES SUGGESTED IN REFERENCE 8 
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TABLE C-II. - NONSYSTEMATIC CODES SUGGESTED IN REFERENCE 9 




I 

o 


ct 
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Originally proposed by J. A. Heller (ref. 
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(c) Code 3 (V = 2, K = 5) . (d) Code 4 (V = 2, K - 6) . 




(i) Code 9 (V = 3, K = 5) . (j) Code 10 (V - 3, K = 6) . 

Figure C-l. - Shift register representation of the encoders using codes suggested 

in reference 8. 
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(e) Code 26 (V = 3, K = 7). 


(f) Code 27 (V = 3, K = 8) . 


Figure C-3. - Shift register representation of the encoders using codes suggested 

in reference 10. 





APPENDIX D 

SIMULATION DATA FOR VITERBl DECODER SIMULATIONS 


The input and output data for the Univae 1108 computer simulations of the Viterbi 
decoding algorithm are summarized in this appendix. For each code investigated, 
tables D-I to D-III contain the following information. 


Code 

Q 

K 

V 

Information bits 
Output errors 


Identification (the particular code) 

Receive quantization levels (Q = 2 is 2-level or hard-decision; 
Q = 8 is 8-level or 3-bit soft decision.) 

Constraint length of the code 

Number of modulo 2 adders in the encoder (Code rate = 1/V.) 

Number of information bits processed (The number of channel 
bits is V times this number.) 

Number of information bits in error after decoding 


Input error 
probability 


Channel bit error probability at decoder input or information bit 
error probability before coding 


Output error 
probability 


Information bit error probability at decoder output 

/ number of output errors \ 

V number of information bits / 


V N 0 


Ratio of energy per information bit to single-sided noise spectral 
density at decoder input 
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TABLE D-I. - VITERBI HARD-DECISION (Q = 2) SIMULATION DATA 


(NONSYSTEMATIC CODES) a 


Code 

number 

K 

V 

Information 

bits 

Output 

errors 

Input error 
probability 
(b) 

Output error 
probability 

V N 0’ 

dB 

3 

5 

2 

66744 

22 

3.038 xlO" 2 

3. 3 x 10' 4 

5.99 





92 

4. 523 x 10" 2 

1. 38 x 10“ 3 

4.58 





128 

4.981 x 10' 2 

1.92 xlO -3 

4. 32 





915 

6.975 x 10‘ 2 

1. 371 x 10' 2 

3. 38 

9 

5 

3 

66744 

0 

5.00 x 10' 2 

_ _ 

6.09 





11 

6. 993 x 10 -2 

1. 65 x 10 -4 

5. 14 





52 

9. 00 x 10' 2 

7. 8 x 10‘ 4 

4. 31 





127 

10.008 x 10" 2 

1.9 x 10‘ 3 

3.92 

16 

8 

2 

66744 

22 

4. 523 x 10" 2 

3. 3 x 10' 4 

4. 58 
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4. 881 x 10 -2 

6.9 x 10" 3 

4. 32 





628 

6.975 x 10“ 2 

9. 41 x 10' 3 

3. 38 





2839 

8.99 x 10' 2 

4. 254 x 10‘ 2 

2. 55 

21 

8 

3 

66744 

3 

8.013 x 10" 3 

4. 5 x 10" 5 

4. 72 





22 

9.0 x 10" 2 

3. 3 x 10~ 4 

4. 31 





63 

10.008 xlO" 2 

9. 4 x 10 -4 

3.92 





129 

10. 996 x 10" 2 

1.93 x 10" 3 

3. 52 


3. 

Refer to figure 11 and table I. 

Channel bit error probability (after coding) . 
















TABLE D-H. - VITERBI SOFT-DECISION (Q = 8) SIMULATION DATA (NONSYSTEMATIC CODES) 3 


Code 

K 

— 

Information 

— 

Output 

Input error 

Output error 

V N 0’ dB 

number 

H 

bits 

errors 

probability 

probability 

2 

4 

2 

66744 

4 

b 1.0 x 10' 2 

5.99 x 10' 5 

4. 325 





128 

b 2.0 

1.91 x 10' 3 

3.24 





446 

b 3. 0 

6.68X10' 3 

2. 48 

3 

5 

2 

66744 

1 

b 1.0 

1. 49 x 10' 5 

4. 325 





100 

b 2.0 

1. 5 x 10' 3 

3.24 





306 

b 2. 5 

4. 58 x 10' 3 

2.84 





319 

b 3. 0 

4. 78 x f0' 3 

2. 48 

4 

6 

2 

66744 

56 

b 2.0 

8. 39 x io' 4 

3.24 





154 

b 2. 5 

2. 307 x 10' 3 

2. 84 





305 

b 3. 0 

4. 57 x 10* 3 

2. 48 

5 


2 

66744 

16 

b 2. 0 

2.397X 10' 4 

3.24 


mi 



67 

b 2. 5 

1. 004 x IO* 3 

2.84 





190 

b 3. 0 

2.847 X 10' 3 

2. 48 

16 

H 

2 

66744 

8 

b 2.0 

1.2 x 10' 4 

3.24 





56 

b 2. 5 

8. 39 x 10' 4 

2.84 


■ 



167 

b 3. 0 

2. 5 x io' 3 

2.48 

8 


3 

66744 

4 

b 1.0 

5.99 x io' 5 

4. 325 





67 

b 2.0 

1. 004 x 10 3 

3. 24 





224 

b 3.0 

3. 356 x 10' 3 

2. 48 

9 

5 

3 

66744 

4 

C 11 

5.99 x 10' 5 

3. 54 





18 

b 2. 0 

2 . 7 x io' 4 

3.24 





71 

b 2.5 

1. 06 x 10' 3 

2. 84 





189 

b 3. 0 

2. 83 x io' 3 

2. 48 

19 

6 

3 

66744 

22 

b 2. 0 

3. 296 x 10' 4 

3. 24 





55 

b 2.5 

8.24X10' 4 

2.84 





148 

b 3. 0 

2. 217 x 10' 3 

2. 48 

20 

7 

3 

66744 

15 

b 2.0 

2.247X 10' 4 

3. 24 





23 

b 2.5 

3. 446 x 10' 4 

2.84 





65 

b 3. 0 

9. 74 x io' 4 

2. 48 

21 

8 

3 

66744 

15 

b 2. 0 

2. 2 x 10' 4 

2. 72 





64 

b 3. 0 

9. 6 x 10‘ 4 

2. 48 


a Refer to figures 12 to 14 and tables I and II. 
^Information bit error probability (before coding) . 
C Channel bit error probability (after coding). 
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TABLE D-in. - VITERBI SOFT-DECISION (Q - 8) SIMULATION DATA 

(SYSTEMATIC CODES) a 


Code 

number 

K 

V 

Information 

bits 

Output 

errors 

Input error 
probability 
(b) 

Output error 
probability 

V N o- dB 

22 

5 

2 

66744 

2 

3. 0 x 10 -2 

3 X10" 5 

5. 49 





13 

4. 8 

1.9 x 10" 4 

4.86 





43 

5.0 

6. 4x 10" 4 

4. 33 





177 

7.0 

2.65 x 10 -3 

3. 38 





437 

9.0 

7.09 x 10" 3 

2. 55 

25 

5 

3 

66744 

6 

7.0 

9 x 10" 5 

5. 14 





21 

9.0 

3. 1 x 10' 4 

4. 31 





72 

11.0 

1.08 x 10' 3 

3. 54 


a Refer to figures 13 and 14. 
^Channel bit error probability. 
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